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1 Introduction 

Let P be a bounded domain (i.e., a bounded connected open set) in of class N >2. 
We consider the following Steklov problem for the biharmonic operator 

— rAlt = 0, in P, 

0 = 0, ondP, (1.1) 

- ^ = \u, on 5P, 

in the unknowns u (the eigenfunction), A (the eigenvalue), where r > 0 is a fixed positive 
constant, v denotes the outer unit normal to 5P, div^n denotes the tangential divergence 
operator and D^u the Hessian matrix of it. For N = 2, this problem is related to the study of 
the vibrations of a thin elastic plate with a free frame and mass concentrated at the boundary. 
The spectrum consists of a diverging sequence of eigenvalues of finite multiplicity 

0 = Ai(P) < A2(P) < • • • < Aj(P) < • • • , 

where we agree to repeat the eigenvalues according to their multiplicity. We note that problem 
dm is the analogue for the biharmonic operator of the classical Steklov problem for the Laplace 
operator, namely 

f All = 0, in P, 

\ = All, on dP, ^ 

which models the vibrations of a free membrane with mass concentrated at the boundary. 
Problem (EU was first considered by Steklov in [39], where the author provided a physical 
derivation (see also [34]). We refer to [30] for related problems, and to [20] for a recent survey 
on the subject. 

In this paper we are interested in the dependence of the eigenvalues Aj(P) of problem (II.1|) 
on the domain P. Domain perturbation problems have been widely studied in the case of the 
Laplace operator subject to different homogeneous boundary conditions (Dirichlet, Neumann, 
Steklov, etc.), in particular for shape optimization problems. We recall for instance the cele¬ 
brated Faber-Krahn inequality, which says that the ball minimizes the first eigenvalue of the 
Dirichlet Laplacian among all domains with fixed measure (see [111129]). Similar results have 
been shown also for other boundary conditions (see e.g., [TUllUllll]). As for the biharmonic 
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operator, much less is known. Lord Rayleigh conjectured that the ball minimizes the funda¬ 
mental tone of the clamped plate (i.e., the first positive eigenvalue of the biharmonic operator 
with Dirichlet boundary conditions) among open sets with the same measure. This has been 
proved by Nadirashvili m for N = 2, and soon generalized by Ashbaugh and Benguria [6] for 
N = 3, while the general case remains an open problem (see also [551110] )■ Regarding Neumann 
boundary conditions, Chasman m proved that the ball is a maximizer for the fundamental 
tone. We refer to m for a general approach to domain perturbation problems (see also |25|L 
and to |26j for a comprehensive discussion on eigenvalue shape optimization problems for elliptic 
operators. We also refer to [HI [13] where the authors prove analyticity properties in the spirit 
of |31] for Dirichlet and intermediate boundary conditions respectively, and show that balls are 
critical domains for all the elementary symmetric functions of the eigenvalues. 

Problem dni) should not be confused with other Steklov-type problems already discussed 
in the literature. For example, in m the authors consider the following problem 

{ A^u = 0, in D, 
u = 0, on dil, 

Au = X^, oiidn, 

which has a rather different nature. We note that, broadly speaking, one may refer to Steklov- 
type boundary conditions for those problems where a spectral parameter enters the boundary 
conditions. 

The aim of the present paper is to discuss the Steklov problem dm as the natural fourth 
order version of problem (O. We derive problem (EB) starting from a physical model and 
study the relationship with the Neumann eigenvalue problem for the biharmonic operator con¬ 
sidered in [M] . Then we adapt the arguments used in [ST] [32] in order to show real analyticity 
of the elementary symmetric functions of the eigenvalues of dO) and compute Hadamard-type 
formulas, which are used to prove that balls are critical domains. For completeness, we do the 
same also for the Neumann problem as stated in M- Finally, we study problem dm when D 
is a ball and identify the fundamental tone and the corresponding modes (the eigenfunctions). 
By following a scheme similar to that used in [55] we prove that the ball is a maximizer for the 
first positive eigenvalue of problem dni) among all bounded domains of class C^. 

The paper is organized as follows. In Section [5] we derive problem (11.11) providing a physical 
interpretation. In Section [5] we characterize the spectrum and show that problem (11.11) is 
strictly related to the Neumann eigenvalue problem as described in [14]. As a bypass product, 
we provide a further phyisical justification of dnD- In Section [4] we compute Hadamard-type 
formulas and prove that balls are critical domains under measure constraint for the elementary 
symmetric functions of the eigenvalues of dm and of the corresponding Neumann problem 
()3.10l) . In Section [S] we prove the isoperimetric inequality for the fundamental tone. Finally, in 
Section [5] we provide some remarks on problems (11.11) and (13.101) when r = 0. 


2 Formulating the problem 

In this section we provide a physical interpretation of problem dm for TV = 2, which arises in 
the theory of linear elasticity, in particular in the study of transverse vibrations of a thin plate. 
Actually, in Sections [2] and [3] we will consider a slightly more general version of problem (11.11) , 
namely 

{ A^u — tAu = 0, in D, 

0=0, on on, (2.1) 

r- divan {D^u.iy) - ^ = Xpu, on SD, 

where a positive weight p G L°°(9r2) appears in the boundary conditions. The weight p has 
the meaning of a mass density. We shall always assume that r is a fixed positive real number. 
We recall that the tangential divergence divan.F of a vector field F is defined as divanF’ = 
divF|gjj — {DF.v) ■ z/, where DF is the Jacobiam matrix of F. 

As usual, we assume that the mass is displaced in the middle plane of the plate parallel to 
its faces. When the body is at its equilibrium it covers a planar domain D with boundary dVL 
in We describe the vertical deviation from the equilibrium during the vibration of each 
point (x,y) of at time t by means of a function v = v{x,y,t). We suppose that the whole 
mass of the plate is concentrated at the boundary with a density which we denote by p{x,y). 
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Moreover, we assume that p{x, y) is bounded and positive on dfl. Under these assumptions, 
the total kinetic energy of the plate is given by 

\ f Pv'^da, 

^ Jon 

where we denote by v the derivative of v with respect to the time t, and by da the surface 
measure on dfl. Now we obtain an expression for the potential energy of the plate. By following 
|431 §10.8], under the assumption that the strain potential energy at each point depends only 
on the strain configuration at that point and that the Poisson ratio of the material is zero, we 
have that the strain potential energy is given by 

W + vly + 2vly)dxdy. 

^ JQ. 

Besides 14, we have another term of the potential energy due to the lateral tension 

= ^ / {vl+vl)dxdy, 

^ Jn 

where t > 0 is the ratio of lateral tension due to flexural rigidity. The Hamilton’s integral of 
the system is given by 



According to Hamilton’s Variational Principle, the actual motion of the system minimizes such 
integral. Let v{x,y,t) be a minimizer for TL. By differentiating (12.21) it follows that v satisfies 


•C2 p p<^2 p 

/ rjpvdadt — p (A^z; — tAv) dxdydt 

1 J oo Jti Jq 

f dv / 2 \ dAv\ 


f dp d'^v 
Iqq dv dv'^ 


for all p G C^(fl X [ti,t 2 ]) such that p{x,y,ti) = p{x,y,t 2 ) = 0 and {x,y) G fl. Since p is 
arbitrary we obtain 


A'^v — tAv = 0, in H, 

0 = 0, on au, 

(D'^v.v) - ^ = 0, on dil 


(2.3) 


for all t G [^ 1 ,^ 2 ]- We remark that we wrote div^n (D'^v.v) instead of divgo qq since 

(^D'^v.v)g^ = D'^v.v — §^v and ^ = 0 on dil. 

We separate the variables and, as is customary, we look for solutions to problem (12.31) of the 
form v{x,y,t) = u{x,y)w{t). We hnd that the temporal component w(t) solves the ordinary 
differential equation —w{t) = Xw{t) for all t G [^ 1 ,^ 2 ], while the spatial component u solves 
problem (EH). 


3 Characterization of the spectrum. Alternative deriva¬ 
tion of the problem 

In this section we prove that the spectrum of the eigenvalue problem ()2.1|) is discrete. In 
particular, each eigenvalue is non-negative and has finite multiplicity and there exists a Hilbert 
basis of the standard Sobolev space of eigenvectors. Then we provide a further derivation 

of problem (12.11) . Namely, we show that this problem can be seen as a limit of eigenvalue 
problems for the biharmonic operator with Neumann boundary conditions and mass density ps 
which concentrates in a neighborhood of the boundary as e goes to zero. We refer to mm 
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for similar discussions concerning second order problems. We observe that the asymptotic 
analysis of mass concentration problems for second order operators has been performed by 
several authors by exploiting asymptotic expansions methods, see e.g., [331 [53] a-nd the references 
therein. We also mention the alternative approach based on potential theory and functional 
analysis proposed in [HISS]. 

Note that here and in the sequel we shall not put any restriction on the space dimension. 
Thus fl will always denote a bounded domain in of class C^, with N >2. 


3.1 Analysis of the spectrum of problem (12.Ih 

Let p G TZ^, where TZ^ := {p G L°°{dn) : essinfa;g 9 n/o(a;) > 0}. We consider the weak formula¬ 
tion of problem in, 


/ D'^u:D'^(p + T'Vu-'V(pdx = X pupda, \/(p € , (3.1) 

Jn Jan 


in the unknowns u G A G M, where 


N 


D^u : D'^p = 


i,i=l 


d'^u d'^p 
dxidxj dxidxj 


denotes the Frobenius product. Actually, we will obtain a problem in since we need to 

get rid of the constants, which generate the eigenspace corresponding to the eigenvalue A = 0. 
We denote by Jp the continuous embedding of into H^{V,y defined by 




We set 


[u][(/9] := f pupda, Vu G L^{dV,),p G 

Jan 

Hp’^{Q) := ju G : f puda = ol , 

I Jan j 


and we consider in the bilinear form 


<u,v>= / D^u : D'^v + tVu ■ Wvdx. (3-2) 

Jn 

By the Poincare-Wirtinger Inequality, it turns out that this bilinear form is indeed a scalar 
product on Hp’^{^l,) whose induced norm is equivalent to the standard one. In the sequel 
we will think of the space iLp’‘^(n) as endowed with the form (13.2|) . Let F{Q) be defined by 
F{id) := {G G iL^(r2)' : G[l] = O}. Then, we consider the operator Vp as an operator from 
F[p''^{n) to F{D,), defined by 



D^u : D^p + tVu ■ Vpdx, Vit G p G F[^{VL). 


(3.3) 


It turns out that Vp is a homeomorphism of onto F{^1). We define the operator tt^ 

from ij2(II) to by 



fan P^dcr 
fan ' 


(3.4) 


We consider the space F[^ (SI)/R endowed with the bilinear form induced by (13.21) . Such bilinear 
form renders (SI)/R a Hilbert space. We denote by the map from Fl^ (SI)/R onto (SI) 
defined by the equality tt^ = o p, where p is the canonical projection of F[^{^) onto 
iJ^(SI)/R. The map turns out to be a homeomorphism. Finally, we define the operator 
Tp acting on iL^(SI)/R as follows 


Tf := (4’'")-' o (Pf )-i o o W o 


where Tr denotes the trace operator acting from iJ^(SI) to L^(cISI). 


(3.5) 
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Remark 3.6. We observe that the pair {X,u) of the set (R \ {0}) x \ {0}) satisfies 

iS.l\) if and only if \ > 0 and the pair (A“^,p[u]) of the set M x ((R^(ri)/]R) \ {0}) satisfies the 
equation 

X~^p[u] = T^p[u]. 

We have the following 

Theorem 3.7. The operator is a non-negative compact selfadjoint operator in R^(n)/]R, 
whose eigenvalues coincide with the reciprocals of the positive eigenvalues of problem SS.1\) . In 
particular, the set of eigenvalues of problem LI. 1\) is contained in [0,+oo[ and consists of the 
image of a sequence increasing to +oo. Each eigenvalue has finite multiplicity. 

Proof. For the selfadjointness, it suffices to observe that 

< >// 2 (n)/R = < o o oTro7r;“’‘^M,r; >^ 2 (Q)/R 

= 'Pp ° Jp o Tr o [J/v] 

= f7’^[Tr o 7rJ’‘^'u][7r*’‘^?;], Vu, z; G i7^(n)/R. 

For the compactness, just observe that the trace operator Tr acting from to Lfi{dLl) is 

compact. The remaining statements are straightforward. □ 

As a consequence we have that the spectrum of is of the form 


0 < Ai < A 2 < • • • < Aj < • • • 


Note that the first positive eigenvalue is A 2 as proved by the following 

Theorem 3.8. The first eigenvalue Ai of i3.1\} is zero and the corresponding eigenfunctions 
are the constants. Moreover, X 2 > 0. 


Proof. It is straightforward to see that constant functions are eigenfunctions of dsn with 
eigenvalue A = 0. Suppose now that u is an eigenfunction corresponding to the eigenvalue 
A = 0. Then we have 


/ + rlVupda; = 0, 

Jn 


where i ox dx )^■ Sin.ce Vit = 0, it follows that u is constant. Then the 

eigenvalue A = 0 has multiplicity one. □ 


Thus A 2 is the first positive eigenvalue of (13.11) which is usually called the fundamental tone. 
Note that we can charactrize A 2 by means of the Rayleigh principle 


L \D^u\^ + rlVupdx 
min - 

fan 


fan 


(3.9) 


3.2 Asymptotic behavior of Neumann eigenvalues 

We consider the following eigenvalue problem for the biharmonic operator subject to Neumann 
boundary conditions 

{ A^m — tAu = Xpu, in f2, 

0=0, on an, (3.10) 

- ^ = 0, on an, 

where p G := {p G L°°{Tt) : ess infa;gop(a:) > 0} (we refer to [T3] for the derivation of the 
boundary conditions). It is well known that this problem arises in the study of a free vibrating 
plate whose mass is displaced on the whole of n with density p. 

Let us denote by n^ the set defined by 

n^ := {x G n : dist(a;, an) > e}. 

We fix a positive number M > 0 and choose the family of densities pe defined as follows 
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(3.11) 


Pe{x) 


e, 

M-£|ne| 

|n\ne| ’ 


if X € He, 
if a: € H \ He, 


for £ G]0,eo[, where eo > 0 is sufficiently small. If in addition we assume that H is of class C^, 
Eq can be chosen in such a way that the map x ^ x — ve \s & diffeomorphism between i9H and 
9He for all e €]0 ,Eo[. We note that J^p^dx = M for all e €]0,eo[. We refer to the quantity M 
as the total mass of the body. 

We prove, under the additional hypothesis that H is of class C^, convergence of the eigenval¬ 
ues and eigenfunctions of problem (I3.10|) with density pe to the eigenvalues and eigenfunctions 
of problem with constant surface density when the parameter e go to zero (see Corol¬ 
lary [2211) • This provides a further interpretation of problem (13.11) as the equation of a free 
vibrating plate whose mass is concentrated at the boundary in the case of domains of class . 

Problem (13.101) has an increasing sequence of non-negative eigenvalues of finite multiplicity 
and the eigenfunctions form a Hilbert basis of iJ^(H). We consider the weak formulation of 
problem (13.101) with density pe. 


/ D^u ■. D^if + tVu-V ipdx = X / peUtpdx , Wif G H^{^), (3-12) 

Jq Jq 

in the unknowns u G A G K. In the sequel we shall recast this problem in i7^(H)/]R 

since we need to get rid of the constants, which generate the eigenspace corresponding to the 
eigenvalue A = 0. We denote by i the canonical embedding of iJ^(H) into i^(H). We denote 
by the continuous embedding of L^(H) into , defined by 




-AT 


We set 


[u][(p] := / peupdx, Vu G L^(H),(p G 

Jn 

:= G i?^(H) : J up,dx = o| . 


Such form 


In the sequel we will think of the space as endowed with the form 

defines on a scalar product whose induced norm is equivalent to the standard one. 

We denote by the map from iJ^(H) to Hp’^(H) defined by 



In 

In Pedx ’ 


for all u G We denote by 7:^1^ the map from iJ^(H)/R onto {VL) defined by the 

equality tt^ = ttJW op. As in (13.31) . we consider the operator as a map from iJp’-^(H) to 
F{Vl) defined by 

V^^[u][ip\ := / D^u : D^ip -b tVm • V(pdx, Vu G p G 

" Jn 

It turns out that is a linear homeomorphism of onto T'(H). Finally, let the operator 

from i7^(H)/IR to itself be defined by 


■= (tt^ 
Pb \ p 


i'p] 




(3.13) 


Remark 3.14. We observe that the pair {X,u) of the set (R\ {0}) x (i?p’^(H) \ {0}) 
i3.12\} if and only if X > 0 and the pair (A“^,p[rt]) of the set R x ((i7^(H)/R) \ {0}) 
the equation 




satisfies 

satisfies 


As in Theorem 13.71 it is easy to prove the following 


Theorem 3.15. Let LI be a bounded domain in of class and e g]0, Eo[- The operator 
is a compact self adjoint operator in i7^(H)/R and its eigenvalues coincide with the reciprocals of 
the positive eigenvalues Xj{pe) of problem \3.1S\} for all j G N. Moreover, the set of eigenvalues 
of problem LS.12\) is contained in [0,-|-oo[ and consists of the image of a sequence increasing to 
-boo. Each eigenvalue has finite multiplicity. 
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We have the following theorem on the spectrum of problem (j3.12|) (see also Theorem I3.8|l . 


Theorem 3.16. The first eigenvalue Ai(p£) of i 3.1 21) is zero and the corresponding eigenfunc¬ 
tions are the constants. Moreover, \ 2 {,Pe) > 0. 

Now we highlight the relations between problems (13.11) and (13.121) when Q. is of class . In 
particular we plan to prove the following 

Theorem 3.17. Let LI be a hounded domain in M.^ of class . Let the operators T ^m and 

from to itself be defined as in \3. 51) and \3.13\) respectively. Then the sequence 

{^,^}ee]o.eo[ converges in norm to T ^m as e ^ 0. 

We need some preliminary results in order to prove Theorem 13.171 We remark that = 
for all c G K, with c ^ 0. This can be easily deduced by (13.41) . 

Lemma 3.18. Let LI be a bounded domain in of class . Let G TZ'^ be as in i3.11\) . 
Then the following statements hold. 

i) For all ip G H'^{Ll)/'R, —>■ 'k\^[lp\ in L^{Ll) {hence also in H^{Ll)) as e ^ 0; 

ii) If Us ^ u in then (possibly passing to a subsequence) T^^fi^[ue\ 7r5’‘^[it] in 

Lfi{Ll) as e ^ 0; 

Hi) Assume that Ug, u, Wg, w G H^{L1) are such that Ug ^ u, Wg ^ w in L^{L1), Tr[u£] —>■ Tr[rt], 
Tr[r(;e] —^ Tr[ry] in L^{dLl) as e —>■ 0. Moreover assume that there exists a constant C > 0 
such that II Vu£||^ 2 (q) < C, || Vr(;e||^ 2 (Q) < C for all e G]0,£ro[- Then 

/ pe {ug - u) Wgdx 0 

Jn 


and 


as e —>■ 0 



w) udx —>■ 0 


Proof. The proof is standard but long and we omit it. We refer to |34j and references therein 
for details. We remark that in order to prove this lemma we need to assume that Ll is of class 
since the Tubular Neighborhood Theorem is used to perform computations on the strip 
Ll\Llg. □ 

Proof of Theorem \3.17\ It is sufficient to prove that the family of compact oper¬ 

ators, compactly converges to the compact operator T ^m . This implies, in fact, that 

iaf2| 


lim 
€—^0 


( 


rriJ\f 
^ Pe 



= 0 . 


(3.19) 


Then, since the operators and T'm are selfadjoint, property (13.191) is equivalent to 

convergence in norm. We refer to nm] for a more detailed discussion on compact convergence 
of compact operators on Hilbert spaces. We recall that, by definition, compactly converges 
to T^m if the following requirements are fulfilled: 

i) if ||'ae||^ 2 (Q)/R < C for all e G]0,e:o[j then the family {T^Me}£g]o,£o[ relatively compact 

inij2(^l)/M; 

ii) if Me —>■ M in iJ^(H)/R, then T^Ug —)► T^m u in iJ^(f2)/R. 
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We prove i) first. For a fixed u € if^(r2)/K. we have 


lim f p^TTj;^[u]dx = lim f — 7'‘5’‘^[u]') dx 


M 

]MI 


'an 


7rJ’‘^[M]c?(T. (3.20) 


By Lemma [3.181 we have that the first summand in the right-hand side of p.20l) goes to zero as 
e —>■ 0 and by standard calculus we have that the second term goes to zero as e —>■ 0. Moreover, 
the equality o o (Pf )“^ holds. Therefore, it follows that T^u is 

bounded for each u € Thus, by Banach-Steinhaus Theorem, there exists C' such 


that \\T^\ 


< C for all £r e]0,£ro[- Moreover, since ||ue|1 


all e €]0,eo[, possibly passing to a subsequence, we have that ^ u in 
u € H‘^{Vt)/R. This implies that, possibly passing to a subsequence, - 

as e —0, for some w S iL^(fl)/R. We show that w = T*^, 

We ■= T^Ue- By Lemma [3.181 we have 




for some 
w in iJ^(f7)/]R 


M u. To shorten our notation we set 

|sn| 


lim [ K]) : [(/?]) + [lUe]) ■ V(7r“f^ [ip])dx 

= [ + TV(7rf’‘^H) • V(7rf’‘^[</j])da:, 

Jn 

for all ip G H^(fi)/R. On the other hand, since o Wg = {>7^ o * o we have 

that 


K]) : D'^{T^if[‘P]) + rViirlf [w^]) • V(7r#f^ [(f])dx 


Then, by Lemma 13.181 in) we have 



Yp\dx 


dx 


<w,(p >ff 2 (n)/R= lim <We,>p >h‘^{q)/r= lim PeT^p^ [^elTr^f^ [‘p]dx 

= 1™ ^ Pe - 7r?’‘® [u]) Jpf [tp]dx + 1™ ^ M - ^1 

-f lini / PeTri^[u]-K\;^[p\dx =-^ [ TT\;^[u]TT\;^[ip\d(J =<T%_U,ip>H'i(Q,)IR, 

|c'“l dan 

hence w = T ^m u. In a similar way one can prove that ||we||^ 2 (q)/r —t In fact 


lim IkeII 

e-fO 


2 

ff2(0)/R 


lim 
£ —^0 

//'( 

[Ue] - TtPM 

) T^lf[We]dx 



-|- lim 

s —^0 ^ 



^ dx 


-1- lim 

f P£7r}’‘^[u] 

[we] - 7r*’‘^[u;]^ 

1 dx 


£—^0 

Jn ^ 




-I- lim / p£7r?’‘^[M]7r?’‘^[w]dx 


This proves i). As for point ii), let Ug —>■ u in iJ^(0)/R. Then there exists C” such that 
Il'^ellff 2 (f 2 )/R ^ 7!" for all £. Then, by the same argument used for point i), for each sequence 
e, —>■ 0, possibly passing to a subsequence, we have Ue. —>■ T^m u- Since this is true for each 

{e,},gN, we have the convergence for the whole family, i.e., T^Ug —>■ T^m u. This concludes 
the proof. □ 














Now we recall the following well-known result. 

Theorem 3.21. Let H be a real Hilbert space and a family of bounded selfadjoint 

operators converging in norm to the bounded selfadjoint operator A, i.e., lime_>.o ||Ae — 

= 0. Then isolated eigenvalues X of A of finite multiplicity are exactly the limits of eigenvalues 
of Ai;, counting multiplicity; moreover, the corresponding eigenprojections converge in norm. 

Thanks to Theorem 13.211 as an immediate corollary of Theorem 13.171 we have 

Corollary 3.22. Let LI be a bounded domain in of class . Let Xj{pe) be the eigenvalues 
of vroblem on LI for all i € N. Let A.. ?’ N denote the eiaenvalues of vroblem m> 

corresponding to the constant surface density Then lim£_>o Xj{pg) = Xj for all j € N. 

4 Symmetric functions of the eigenvalues. Isovolumetric 
perturbations 

In this section we compute Hadamard-type formulas for both the Steklov and the Neumann 
problems, which will be used to investigate the behavior of the eigenvalues subject to isovol¬ 
umetric perturbations. To do so, we use the so called transplantation method, see m for a 
general introduction to this approach. We will study problems (II.ip and (13.101) in for a 

suitable homeomorphism </>, where LI has to be thought as a fixed bounded domain of class . 
Therefore, we introduce the following class of functions 

<I)(n) = |(^ g (^))^ ■ ^ injective and inf | det D<j)\ > o| . 

We observe that if LI is of class and (j € d)(n), then also (/)(Ll) is of class and g 

$((/)(n)). Therefore, it makes sense to study both problem (11.11) and problem (I3.10p on (/>(L}). 
Moreover, we endow the space C^(L}) with the standard norm 

WfWc^n) = sup jD^^f^x)]. 

|ci:|<2, 

Note that $(0) is open in (C^(n))^ (see [311 Lemma 3.11]). 

We recall that it has been pointed out that balls play a relevant role in the study of isovol¬ 
umetric perturbations of the domain LI for all the eigenvalues of the Dirichlet and Neumann 
Laplacian. We refer to [311133) , where the authors prove that the elementary symmetric func¬ 
tions of the eigenvalues depend real analytically on the domain, providing also Hadamard-type 
formulas for the corresponding derivatives. Then, in [32] they show that balls are critical points 
for such functions under volume constraint. 

From now on we will consider problems (HH) and (13.101) with constant mass density p = 1. 


4.1 The Steklov problem 

We plan to study the Steklov problem in the domain 4>{L}) for (f g $(12), i.e., 

( A^u — tAu = 0, in (f(LI), 

< 0 = 0 , on dm, ( 4 - 1 ) 

I “ div 9 ,p(^n^(B^u.iy) - ^ = Xu, on dcjiLl). 

To do so, we pull it back to LI. Therefore, we are interested in the operator from (LI) 
to F(Ll), defined by 


V^[u][(p] -.= j (D^(uo(f) ^) o (jf) : (D"^ (g} o (j) ^) o c)))\det D(j)\dx 
Jn 

+ T j ex/ (u o (j)~^) o (j)) ■ (V(ip o (j)~^) o (f))\ det D(j)\dx, yu G (LI), ip G H'^(LI), (4.2) 

J Q 


where 


710(12) := g 71^(12) : J m|j^(V0) ^|| det Zl()>|c?cr = o| . 
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Moreover, for every (j) € we consider the map from L^{dfl) to H^{ny defined by 

J?[u][if]-.= f (j))~ ^\\det D(j)\da, yu G L'^{dV,),(f £ 

Jan 

It is easily seen that the form (14.21) is a scalar product on (il). We will think of the space 
as endowed with the scalar product (14.21) . We denote by the map from to 

defined by 

5 Jg^u\v{V \ det D(l)\dx 

/gQ II det ’ 

and by the map from i/^(n)/]R onto {Vl) defined by the equality tt^ = 7r^’‘^op. Clearly, 
TT^’*^ is a homeomorphism, and we can recast problem m as 

X-^u = 

where 

)■' ° ° Tr o Trf. 

The operator Wg can be shown to be compact and selfadjoint, as we have done for the operator 
Tp defined by (j3.5l) in Theorem 13.71 Isee also [33l Theorem 2.1]). 

In order to avoid bifurcation phenomena, which usually occur when dealing with multiple 
eigenvalues, we turn our attention to the elementary symmetric functions of the eigenvalues. 
This is the aim of the following 

Theorem 4.3. Let LI be a bounded domain ofM.^ of class . Let F be a finite non-empty 
subset o/N. Let 


An[F] := {</> G : Xi[fi] i {Xfifi] : j G F} V/ G N \ F}. 

Then the following statements hold. 

i) The set .An[F] is open in 4>(n). The map Pp of Aq[F] to the space C (L.), 

which takes G .4n[F] to the orthogonal projection of Fl^'^{LI) onto its (finite dimensional) 
subspace generated by 

|u G H^’^{Ll) : Pf^[u o = Xj[(j)\J^ o Tr[M o (j)~y for some j G f| 
is real analytic. 

ii) Let s G {1,..., |F|}. The function Ap^s from .4n[F] to defined by 


is real analytic. 

Hi) Let 

0n[F] :={(/) G yIn[F] : \j{4>\ have a common value Vj G F}. 
Then the real analytic functions 



of An[F] to M.^ coincide on 0a[F] with the function which takes (p to Xp[(j)]. 

Proof. The proof can be done adapting that of [33J Theorem 2.2 and Corollary 2.3] (see also 

ED). □ 

In order to compute explicit formulas for the differentials of the functions Ap^s-, we need the 
following technical lemma. 
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Lemma 4.4. Let fl be a bounded domain in of class , and let cj) G be such that 

^{Ll) is of class C^. Let ui,U 2 € IL^(fl) be such that Vi = ut o g for i = 1,2 

and 

d'^vi d‘^V2 ^ 

Then we have 


d\A=A'P^{'^][ui][u2] = [ {D'^vi : D‘^V2 + tVvi ■\7v2)ti'-J^da 

^ Id' n + diva^(n)(-D^W 2 -J^)V?;i) • fida 


f d<^{€l) 

( dl\vi 


f d(j}{VL) 


\ dv 


VV2 


dAv2 

dv 


( dvi 


dv2 


Vvi ■ ^da — T I Vf 2 + -^Vvi • ^da 


' d4>{0.) 


\ dv 


dv 


/ ((A^iii — tAz;i)Vi'2 + {A^V2 — rAw2)Vz;i) • fida, (4.5) 

J 0(0) 


for all Ip g {C'^{n))^, where n = tp o cp 
Proof. We have 


= J id\^^^D^{uiop-^)op)[tp]: {D‘^{u2 0p~^)o4>)\detDp\dx 

+ ^j {d\^^^yiuiop-^)op)['il;]-{V{u2 0p~^)op)\detDp\dx 

+ J {D^{ui o4>-^)op) : {d\^^^D‘^{u2 0p~^)op)[tp]\detDp\dx 

+ tJ {V{ulOp-'^)op)■{d\^^^V{u2 0p~^)op)[^p]\detDp\dx 
+ J {D‘^{ui o o p) : {D‘^{u2 o P~^) o p)d\^^^\del Dp\[%p]dx 

+ '^j {'^iuiop-'^)op) ■ (y{u2 0(p~'^)op)d\^^^\detDp\['ijj]dx, (4.6) 
and we note that the last two summands in (14.61) equals 



[D^vi : D^V 2 + rVfi • Vf 2 ) divfidy. 


(See also Proprosition l4.18l) . By standard calculus we have (see [TSj formula (2.15)]) 


D^{u o p-^) op = {yp)-*D^u{yp)-^ 



du 

dxk 



ij 


where a = (V(/)) This yields the following formula 


N 


dv 


d\^=^{D'^{uop ^)op)[ip]o4> ^ = -D'^vVfi-Vf/D'^v-^—D'^fj.r, (4.7) 

where = ip o (l)~^ and v = u o (j)~^. We rewrite formula (|4.7p componentwise getting 
{d\^^^iD'^iuop-^)op)[p]op-p__ 

N 

= -E 


d'^v dfir d'^v dy.r d'^y.r dv 


dy,dyr dy^ dyjdy^ dyi dyidyj dyr 
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Moreover (see m Lemma 3.26]) 


N 


o r‘) o m] o r'), = -g 

Now we use Einstein notation, dropping all the summation symbols. The first summand of 
the right hand side of (|4.6I) equals 


I^(n) ( 


d^Vi dflr ^ d^Vi dflr 


d^fj,r dvi \ d'^V2 


Uio.) \dyidyr dy^ dyjdyr % dytdyj dyr J dyidyj 

In order to compute (I4.8L integrating by parts, we have 


dy. 


(4.8) 


f d'^vi dyr d‘^V2 
k{Q.) dyidyr dyj dyidyj 


f dvi dyr d^V2 

Jd4,(Q.) 9yi dyj dyidyj 


f dvi ddivy d'^V2 
h{n) dyi dyj dy^dy. 


dy - 


dvi dyr d‘^V2 , 

r'rv ;——da — 


f dvi dyr 

JdMQ) dyi dyj 


f dvi dyr d^V2 
Ln) dyi dyj dyidy^dyr 
f dvi dyr d^V 2 


d4>{Q.) dyi dyj dyidyj 


dy 

dy 


Ia4>{n) dyi 


dvi d‘^V2 , 

awy——;;—Vida 


dyidyj 


l^{Q) dyi dyj dyidyjdyr 
/ D^vi : D^V2divydy 


IcjiQ) 




div/rVui • VAv 2 dy, 


and 


dvi d'^yr d'^V2 


/0(n) dyr dyidyj dyidy^ 


dy= / 


dvi dyr d'^V2 , 

— d(7 


d 4 ,{Q) dyr dyi dyidyj 


d'^vi dyr d‘^V2 
iio.) dyrdyj dyi dyidy^ 
dvi dyr d‘^V2 

-Vi 


■dy - 


da — 


JdHQ) dyr dyi dyidyj 
f dvi dyr d‘^V 2 
la^{n) dyj dyi dyidyj ^ ^ 

dv\ dyr 


+ 

dvi dy. 


1“ dvi dyr dAv2 
Ln) dyr dyi dyi 
f dvi dyr dAv2 
Lf2) dyr dyi dyi 
dvi ddivy d^V 2 
HQ) dyj dyi dyidyj 
d^V2 


dy 

dy 


dy 


0(0) dyj dyi dyidyjdy, 


d\2 , 

—da — 


JdHQ) dyr dyi dyidyj 

f dvi dyr d‘^V 2 

IdHQ) dyj dyidyidyj ^ ^ , 

f dvi d'^V 2 
/ -^div/i Vida - 

Jaipin) dyj dyidyj 


-dy 

f dvi dyr dAv2 
Lf2) dyr dyi dyi 
dvi dyr d^V2 


dy 


0(0) dyj dyi dyidyjdyr 
D'^Vi : D'^V2(\ivydy 


dy 


/0(O) 


/ Vwi • VAu2div/xd?/. 


We also have 


dvi dyr dv2 


dv2 , 


-- —dy = / Vui • yda — 


70(0) dyr dyi dyi 


J d(f){0.) 
dv2 d'^vi 
0(0) dyi dyidyr' 


dv 


/0(O) 


Az; 2 Vui • ydy 


dv2 


yrdy = / ——Vui • yda — 




dv 


/0(O) 


Az; 2 Vui • ydy 


d<p{Q) 


Vui • Vv 2 y ■ vda + 


f dvi d‘^V2 

dyi dyidyr 


yrdy- 


/0(O) 


Vui • 'S/v 2 ddvydy. 
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It follows that 


D^Vi : D^V2dw^dy — 
' dvi 


4(f2) 


dvi d'^V2 dv2 d^vi \ dy. 


JdHQ) V% dyidyj 
d^v2 dv2 d^vi 


dyi dy.dyj 
5/i, 


/^(n) V% dyidyjdyr dyi dy^dy^dy, 
'dvi d‘^V2 dv2 d'^vi 


dy. 


-dy 


Idim ydy, dyidyj ' dyt dy,dyj J 

(Vtii • VAti 2 + Vti 2 • VAwi) divydy 


0(n) 


I 


dvi d'^V2 dv2 d'^vi 


di{,Q.) ydyr dyidyj 
dvi dAv2 


dyr dyidyj 
dv2 dAvi 


dfJ-r , 
Vj-^da 

dy^ 


0(n) \dyr dyi 
(dv, 

d<l){Q) V 


dyr 

dyr dyi J dyi 


dy 


9yj 


rda 


dv 

„ Vv2 + -T^Vvi ) • yda 
ou oiy 


+ r 


L 

/ {AviVV 2 + Av 2 Vvi) ■ ydy + T / Vvi ■Vv 2 y ■ vda. (4.9) 

Ji(n) Jdiin) 


Now we recall that 


dfi ~ 

divy = div 9 ^(f 2 )/i +-^ ■ v on d(l){n), 


(see also [ISl §8.5]) and that, since i/ = V6, where b is the distance from the boundary defined 
in an appropriate tubular neighborhood of the boundary, then = (VvY and = 0, from 
which it follows that 

We will use these identities throughout all the following computations. 

Using the fact that 

i9^ui d'^V2 


du'^ dv'^ 

we get that the sixth summand in (HU) equals 


= 0 on 90(U), 


d(p{Q) 


^ , 9v2 , „2 N 


[dyr^^ dy 


d<p{Q) 


■ V^T 




• da 


dvi d'^V2 , dv2 d'^vi \ dyr 


/ O0(n) 

+/. 


Jdim \9yr dv'^ 

d'^vi d‘^V2 d'^V2 d'^vi 


yrda 


Vjyrda 


lai(Q.)\dyidyr dyidyj ' dyidyr dy^dy^ J ^ 

^ (diva,^(n)(£'^^^i-i^)a,^(a)Vu2 + div9^(j2)(^^^2-*")a0(n)^^i) ’ (4-10) 
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The seventh summand in (j4.9|) equals 

[ f ^Vi;2 + —— [ (A‘^vi'Vv2 + ^^V2'^vi)-^da 
Ja^{Q.) \ ov av J 


'4>(n) 

d'^vi dAv2 d‘^V2 dAvi 


Urdy 


h{n) \dyidyr dyi dy^dyr dyi 

( Vtii^ • pida - [ (A^vi'Vv2 + A^V2'Vvi) ■ y-da 

\ Ol' ov J J^{Q,) 

f 7 f f 9vi d'^Av2 dv2 d'^Avi 

- / {Vvi ■ y Av 2 +yv 2 ■ y Avi) y ■ vda + / I 

J 5(^(0) J V 


d^{Q.) 


yrdy 


0(n) V dyidyr dyi dy^dyr 

+ / (Vt;i • VAW 2 + Vt '2 • VAwi) div/xd?/. (4.11) 
Ji(n) 


The second summand in (j4.9l) equals 


/ V(Vi;i • yV 2 )y{yr)i^rda 

{Vvi-Vv2)V Ql,m{yr)Vrda - [ ^{Vvi ■S/V2)^^VrdcF. (4.12) 

Jaiio.) ov ov 


dcl>(Q.)\ 


J 00 ( 0 ) 

The third summand in (14.91) equals 

' dvi d^V2 dv2 d^vi 


+ 


Vjyrda 


la^in)\9y,dyidyjdyr ' dyi dy.dyjdyr J ^ 

f ( dv\ d'^Av2 dv2 d'^Avi \ 
hi.^) \^yi % dyidyr) ^ 

r / d^vi d^V2 , d‘^V2 d^vi 

Jim V 


+ 


im \dytdyj dyidy^dyr dy^dyj dy^dy^dyr 
dvi d^V2 , dv2 d^vi 

00 ( 0 ) 


yrdy 


him V 

/ 

J 6( 


dyjdyr dy^ dyidyjdyr 
dvi d‘^Av2 dv2 9^Atii 


Vjyrda 


lim'^^yi ^y^^y^ ' ^vi ^vi^yr J 

— / D^vi ■. D^V2y ■ vda + / D^vi : D^V2divydy. (4.13) 

J 00(0) J 0(0) 

From (I4.9F(|4.13L it follows that 

d\,j=^'P^[i’][ui][u2] = ■'^^^2)Va^(fj)(/r^)z/7.dcr 

f d dll f d 

— / T^(Vni • V'i; 2 )-;v^t'rdo- + / —{WVi ■'S/V2)divyda 

Jaim o^ ov Jaim o^ 

d'^vi d'^V2 d‘^V2 d'^vi \ ^ 

dyidyr dy^dyj ^ dyidyr dy^dyj ) 
dvi d^V2 dv2 d^vi \ ^ 

him V% dyidyjdyr dy^ dyidyjdyr) 

— / D^vi : D^V2y ■ vda — / (Vrii • \ 7 Av 2 + Vr;2 • VArii) y ■ vda 

J 00(0) J 00(0) 

+ / ( ^ ^^2 + ^ ■ yda — f (A^vi'Vv2 + A‘^V2'^vi)-ydo' 

Jaim \ ov dv J Jim 


iipm 

f (dvi dv2 ^ \ , 

- T / -^^'*^2 + ■ yda 

Jaim \ov dv 
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+ T / (Az;iVt'2 + Ati2V?;i) •//dy + T / Vvi ■'Vv2fJ. ■ i^da 

J (^(O) J d^{il) 

f 9 

i(n)(Vi'i ■ V?;2)Vg^(fj)(y^)z/^dcr + ’ '^^2)divg^(fj)Mdcr 

^ij'n (^^^a0(O)(^^^i-'")a0(a)'^^2+div9^(f2)(£|2i;2.i/)a^(0)Vi;i) ■ ^ida 


f d^{€l) 


+ 


' d^{Q.) 

( dlS.v\ 


> O0(n) 


V du 


dAv2 ^ , , 

^V2 H-- \/Vi ■ ^da — T 


dv 


d^{0.) 


^Vz;2 + ) • ndcr 

ov ov 


d f d 


-(Vi'i • S/V 2 ) ) firda 


Jaij>{Q.) \dyr 

/ D^vi : D‘^V 2 y ■ vda — / (Viii • VAv 2 + Vv 2 ■ VAvi)y ■ vda 

Jd4>{n) Ja4>{n) 

— / (^A‘^vi'Vv2 + A^V2'^Vi) ■ yda 

J (^(O) 

+ T / {Av\S/V 2 + AV 2 S/V 1 ) ■ ydy + T / Vvi ■'S/v 2 y ■ vda. (4.14) 

J $(Q.) JdiiO,) 


/ (^(n) J 50(f2) 

The first summand on the right hand side of (I4.14|) equals 

/- ^aMn)C^Vl-'^'V2)yvda+ I ^Q^(n){'^'Vl ■^V2) ■ Q^(a)Vr)yrda, 

Jdd>(n) Jd6in) 


' 50(0) 

while the sixth one equals 

/ 50(0) 




/ 50(0) 


■ '^^ 2 ) ■ aHn)’^r-)yrd<7. 


Using the fact that 


Ln, /«<n, 

where K denotes the mean curvature of d(j){il) (see [TOl §8.5]), we obtain 


d\^^^'P^['(lj][ui][u2] = 


50(0) 


^50(0) 


(Vvi • Vn2)/r • vda 


f d /■ 52 

+ / K—{Vvi-Vv2)li-vda+ / —^(Vz;i • Vn2)u • z^dcr 

Jaim ov Jgifm dv^ 


50(0) 


— / : D‘^V2y ■ vda — / (Vr^i • VAv2 + Vn2 • VAvi)fi ■ vda 

J 50(0) J 50(0) 

+ (diva,^(n)(^^^i-'^)a0(n)^^2+ div3^(f2)(i:)2^2.i^)a^(0)V?;i) ■ yda 


+ 


' 50(0) 

/ 9Ar!i 


f 50(0) 


V 9:/ 


dAv2 ^ , , 

Vf2 H—) ■ d'do’ — T 

f 

/0(O) 


' 50(0) 


{dvi.^ 9^2 „ 1 


— / (A^niVv2 + A^n2Vui) • /xdcr 

J 0(0) 

+ T / (Az;iVn 2 + Au 2 V?;i) • ydy + T / \7vi ■ \7v2lJ, ■ vda 

J^{a) Jai{Q.) 

= / A^Vvi ■ VV 2 )y ■ vda — / D^vi : D^V 2 y ■ vda 

Jai(Q.) Jai{Q.) 

— / (Vui • VAn2 + Vn2 • VAni)y • z/du 
v/ 50(0) 
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' d((){Vt) 

( d/^vi 


f 00(0) 


\ dv 


Vt>2 


dAv2 

dv 


{ dvi 


dvo 


Vfi • iida — T I -T^yv2 + -^Vvi • ^da 


00 ( 0 ) 


\ dv 


dv 


/ ((A^iii — tAwi)Vw 2 + (A^ii2 — tAz;2)Vi'i) • ^da 

J 0(0) 


+ r 


' 00 ( 0 ) 


Vz;i • 'Vv 2 fJ- ■ vda. 


Using the equality 

A(Vui • yv2) = VAui • Vv2 + Viii • VAti2 + 2D^vi : D^V2 
we finally get formula (1431) . □ 

Now we can compute Hadamard-type formulas for the eigenvalues of problem mi- 

Theorem 4.15. Let LI be a bounded domain in of class . Let F be a finite non-empty 
subset o/N. Let (p € 0n[U] be such that d(j){L}) G . Let vi^... jriiFl be an orthonormal basis 
of the eigenspace associated with the eigenvalue AF[</i’] of problem |7UP L"^(d(p(Ll)). Then 

\F\ 

_ ! \ H \ _ I \ 

d\^^^iAF,s)[4’] = ~^F W\ 


| U |-1 

s — 


1 E / - 


+ -T|Vwip - ■ vda, 

for all Ip € (C^(U))^, where p = ip o (p^~^\ and K denotes the mean curvature of dp{Ll). 

Proof. First of all we note that vi, ...,r’|F| € H‘^{(p{fl)) (see e.g., [HI §2.5]). We set ui = vi o p 
for Z = 1,..., |F|. For |F| > 1 (case \F\ = 1 is similar), s < |F|, we have 


d‘\rf,=4,i^F,s)['‘P] — ^F M 


| U |-1 
s — 1 


l^’l 


E^l [P{ui)]- (4.16) 


1=1 


We refer to |3TJ Theorem 3.38] for a proof of formula (I4.16|) . 
By standard calculus in normed spaces we have: 


((V ) ^ o {Vf) o o Tr o [ipMui)] [p{ui)] 


= pG 


Now note that: 


1 o o d|^^^ [j^ O Tr o Mpiui)] 


[Pi^i)] 




d| 0=0 (( 7 ^ 0 '^) ^ o i'Pt) M o o Tr o fY[p{ui)] [p{ui)] . 


V'i 


1 o (v?'^ o d|^^^ (j-/ O Tr o bPMui)] 


[piui)] 


[ (kv^+ p-vda - f V{vf)-pda, 

JdMQ.) V OV J JdMQ.) 


(see also [501 Lemma 3.3]) and 


Vp 


^ o i'Pf) [V'] o J-| o Tr o 7r^’‘^[p(Mz)] [p(u;)] 


= [ip][ui][Ff{ui)]. 

(We refer to [551 Lemma 2.4] for more explicit computations). Using formula (|4.5p we obtain 




^ o {V^) bP] o J7| o Tr o 7r^’‘^[p(wi)] b(wi)] 


f [\D'^vi\'^ + t\\ 7 vi\'^) p ■ vda + ( \ 7 {vf)-pda. 

^50(0) J di(Q) 




This concludes the proof. 


□ 
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Now we turn our attention to extremum problems of the type 

min or max 

V(0)=const. ’ V((/>)=const. 

where V{(t>) denotes the measure of <^(n), i.e., 

V((/>) := j dx = I I det £)(/>|c?x. (4-17) 

^ 0 ( 0 ) Jn 

In particular, all t/i’s realizing the extremum are critical points under measure constraint, i.e., 
Ker(iV((^) C Ker dAi?_s[(^]. We have the following result (see [3H Proposition 2.10]). 

Proposition 4.18. Let Lt be a bounded domain in of class . Then the following state¬ 
ments hold. 

i) The map V from $(11) to R defined in J7| ) is real analytic. Moreover, the differential of 
V at (j) G $(11) is given by the formula 

d\^=^V{(t>)['f]= [__ div{tp o ^-^)dy = [ {tf o ■ vda. 

J0(n) Jd4>(n) 


a) For Vo £]0,+oo[, let 

V(Vo) := {<(. G $(11) : V{4>) = Vq}. 

IfV{Vo) 0, then V(Vo) is a real analytic manifold of (C^(ll))'^ of codimension 1. 

Using Lagrange Multipliers Theorem, it is easy to prove the following 

Theorem 4.19. Let Tl be a bounded domain in of class . Let F be a non-empty finite 
subset o/N. Let Vq G]0, +oo[. Let f) gV (Vq) be such that d(f{Ll) G and Xj [(j)] have a common 
value Af[^] for all j G F and A; [(j)] ^ Ai?[())] for all I G N\F. For s = 1,..., |F|, the function (f 
is a critical point for Ap^s on V (Vq) if and only if there exists an orthonormal basis Vi,... ,v\p\ 
of the eigenspace corresponding to the eigenvalue A p {(j)] of problem eip in L^{d(j){^l)), and a 
constant c G R such that 

- r\yvif - \D^vif'^ = c, a.e. on 9^(11). (4.20) 


Now that we have a characterization for the criticality of (j), we may wonder whether balls 
are critical domains. This is the aim of the following 

Theorem 4.21. Let H be a bounded domain o/R^ of class . Let (j) G $(11) be such that 
0(11) is a ball. Let A be an eigenvalue of problem in 0(11), and let F be the set o/j G N 
such that Aj[0] = A. Then Ap^a has a critical point at 0 on U(V(0)), for all s = 1,..., |F|. 


Proof. Using Lemma [4.221 below and the fact that the mean curvature is constant for a ball, 
condition (|4.20l) is immediately seen to be satisfied. □ 


Lemma 4.22. Let B he the unit ball in R^ centered at zero, and let A be an eigenvalue of 
problem ) in B. Let F he the subset o/N of all indeces j such that the j-th eigenvalue of 
problem If.l ) in B coincides with A. Let vi,... ,r|F| he an orthonormal basis of the eigenspace 
associated with the eigenvalue A, where the orthonormality is taken with respect to the scalar 
product in L^{dB). Then 

|F| \F\ |F| 


i=i 


i=i i=i 


are radial functions. 


Proof. Let On{^) denote the group of orthogonal linear transformations in R'^. Since the 
Laplace operator is invariant under rotations, then Vk o A, where A G On(]&.), is still an eigen¬ 
function with eigenvalue A; moreover, {vj o A : j = 1,..., |F|} is another orthonormal basis for 
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the eigenspace associated with to A. Since both {vj : j = 1,..., |J^|} and {vjoA : j = 1,..., |F|} 
are orthonormal bases, then there exists i?[^] G OAr(]R) with matrix (i?ij [A])^ j=i such that 

l^’l 

1=1 


This implies that 



i=i 


li"! 


i=i 


from which we get that i® radial. Moreover, using standard calculus, we get 


|F| |F| |F| 

= E Rjhmi 2 {A]iVvi,oA) ■ (Vvi^oA) = O A|2, 

j=l il.i2 = l 1=1 

and 


D^Vj ■ D^Vj 


\F\ 

E Rjh [A]Rji, [A] A* • [D^i, O A) • A • A‘ • {D\i, o A) • A 

h ,^2 — 1 


= E [^]Rjh [A]A* ■ (D^Vh o ■ (D^Vh o A) • A, 

,^2 — 1 


hence 


\D^v, 


\F\ 

■ D^Vj) = E Rjh [A]Rji2 [A](T'^w;i o A) : o A), 

h ,^2 — 1 


from which we get 


l^’l \F\ 


i=i 


3 = 1 


□ 


4.2 The Neumann problem 

As we have done for the Steklov problem, we study the Neumann problem in i-6-i 

{ A^it — tAu = Xu, in ^(O), 

0 = 0, on dcjiirt), (4.23) 

- diva^(n)(£>^M.F) - ^ = 0, ona(/)(n). 

We consider the operator from (Q) to F{fl), defined by 


[u][(p\ -.= f {D^(uo(j) o (j)) : [D'^ [tp o (j) o ^)| det £>(/)|da: 

Jo. 

+ T I (S/{u o o (j)) ■ {y{(p o o (j))\det D(j)\dx, Vm G iJ?’^(r2), (/? G (4.24) 

Jn 


where 


i?^’^(n) := G Ar^(f2) : J u| det I?())|dx = oj’, 

Moreover, for every 4> G $(f2), we consider the map from L^(n) to H^(D,y defined by 

jyf\u\[ip\:= f u(f \ det D (j)\da, yu G L^{V,),(p G 

Jn 
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We will think of the space {VL) as endowed with the scalar product induced by (j4.24[) . We 
denote by the map from H'^{Vl) to [Vl) dehned by 


TT^ iu) =:u- 


m| det D(/)ldx 
/q I det D(j)\dx 


and by tt^’^ the map from H‘^{n)/M. onto (fl) dehned by the equality op. 

Clearly, is a homeomorphism, and we can recast problem ()4.23|) as 


X-^u = Wfu, 

where := o (V^)~^oJ^ ozott^’^ and i is the canonical embedding of into 

Lp'{Vl). An analogue of Theorem 031 can be stated also in this case. Therefore, we can compute 
Hadamard-type formulas for the Neumann eigenvalues. This is contained in the following 

Theorem 4.25. Let LI be a bounded domain in of class . Let F be a finite non-empty 
subset ofN. Let (p € ©nlf] be such that d(j){L}) € Let vi,... be an orthonormal basis 
of the eigenspace assoeiated with the eigenvalue of problem /14-^dl) in L'^{(p{Ll)). Then 


\p\ 

d\^^^{^F,s)bP] = XI / {XFvf-T\\/vi\'^-\D^vif)p-iyda, 

JafiQ) 


for all Ip £ ((7^(0))^, where p = ip o (p 


Proof. The proof is similar to that of Theorem 005] 

First of all we note that, by elliptic regularity theory, vi,...,v\p\ £ H'^{p{Ll)) (see [TU 
§2.5]). We set u; = ti; o for / = 1,..., |F|. For |F| > 1 (case |F| = 1 is similar), s < |F|, we 
have 

l-^l 

du=^{^F,sm = \dU=pKi^Mui)]\ [pM] . 


1=1 


By standard calculus in normed spaces we have: 


V- 


Af 


dU=^((4’^)-'o«) \j- 


_ -pF 


■M ■ 
oio 7r»’ 

-1 


) Wi[p{ui)\ [p{ui)\ 




[piui)] 


Now note that 


V- 


N 


V- 






d|0=0 ((4’^) ^ o 0*0 T^Y[p{ui)\ b(**;)] • 


[Piui)] 


0 (n) 


vfdivpdy, 


(see also Proposition 14. 181) and 


V- 


AT 


^ o (V^) [pj]oJ^ oio Tr^J^[p{ui)] [p{ui)] 


= OTT^) [V’jbzjlTT^(**;)]■ 


Using formula (14.5|) we obtain 


r- 


AT 


^ o (V^) [ip]oJ^ oio Fh-^[p{ui)] [p{ui)] 


/_ {\D'^vi\'^ + t\\7vi\'^) p-Fda + f_ y{vf)- 

J 50 (f 2 ) J (^(n) 


pdy. 
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To conclude, just observe that 



V(uf) ■ ^^dy= [ 
Ja 


d4>{Q) 


vf /i • vda 



{vf)divfidy. 


□ 


Now we can state the analogue of Theorem 14.191 for problem ()4.23p . 

Theorem 4.26. Let LI be a bounded domain in of class . Let F be a non-empty finite 
subset ofN. Let Vo G]0, +oo[. Let 0 G V(Vo) be such that d(j){Ll) G and Xj\(f] have a common 
value for all j G F and Az[0] ^ Ai;’[0] for all I G N\F. For s = the function 0 

is a critical point for Ap^s on V (Vq) if and only if there exists an orthonormal basis fi,..., 
of the eigenspace corresponding to the eigenvalue Ai;’[0] of problem i4.23\ ) in L^(0(f2)), and a 
constant c G R such that 

{Xpvf — T|Vwip — \D‘^vi\'^) = c, a.e. on 00(17). 

1=1 

We observe that Lemma 14.221 holds for problem (14.2311 as well, since in the proof we have 
only used the rotation invariance of the Laplace operator. Then, we are led to the following 

Theorem 4.27. Let LI be a bounded domain in R^ of class C^. Let 0 G 4*(17) be such that 
0(17) is a ball. Let X be an eigenvalue of problem 114 ■ 231) in 0(17), and let F be the set of j G N 
such that Aj[0] = A. Then Ap^g has a critical point at 0 on V(V(0)), for alZ s = 1,..., \F\. 


5 The fundamental tone of the ball. The isoperimetric 
inequality 

In the previous section we have shown that the ball is a critical point for all the elementary 
symmetric functions of the eigenvalues of problem ini. In this section we prove that the ball 
is actually a maximizer for the fundamental tone, that is 

A2(17) < A2(17*), (5.1) 

where 17* is a ball such that |17| = |17*|. 

5.1 Eigenvalues and eigenfunctions on the ball 

We characterize the eigenvalues and the eigenfunctions of (HID when 17 = i? is the unit ball 
in R^ centered at the origin. It is convenient to use spherical coordinates {r,0), where 6 = 
(01,..., 0jv_i). The corresponding trasformation of coordinates is 


Xi 

= rcos(0i). 


X2 

= r sin(0i) cos(02), 


XN-1 

= r sin(0i) sin(02) • ■ 

' •sin(0Ar_2) cos(0Ar_i) 

Xn 

= r sin(0i) sin(02) • ■ 

' •sin(0jv-2)sin(0jv-i) 


with 01, ...,0vr_2 G [0,7r], 0w-i G [0, 27r[ (here it is understood that 0i G [0,27r[ if iV = 2). 
The boundary conditions of (HD in this case are written as 


d^u 
dr'^ P 


= 0 , 


-fi- =Alt| 

or r-^ ^ \or r j or |r=i 


where As is the angular part of the Laplacian. It is well known that the eigenfunctions can 
be written as a product of a radial part and an angular part (see [14] for details). The radial 
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part is given in terms of ultraspherical modified Bessel functions and powertype functions. The 
ultraspherical modified Bessel functions ii{z) and ki{z) are defined as follows 


ii{z):=z^ In_j^+i{z), 

ki{z) := z^-^K n_^^i{z), 

for I G N, where Iu{z) and K,j{z) are the modified Bessel functions of first and second kind 
respectively. We recall that ii{z) and all its derivatives are positive on ]0,+oo[ (see [U §9.6]). 
We recall that the Bessel functions Jv and Ni, solve the Bessel equation 

z^y"{z) + zy'{z) + - v'^)y{z) = 0, 

while the modified Bessel functions Ii, and solve the modified Bessel equation 

z^y"(z) + zy'{z) + (z^ + v‘^)y{z) = 0. 


We have the following 

Theorem 5.2. Let Ll he the unit hall in eentered at the origin. Any eigenfunetion ui of 
problem E2P is of the form ui{r,9) = Ri{r)Yi{0) where Yi{0) is a spherical harmonic of some 
order I G N and 

Ri{r) = Ay + Biii{y/Tr), 
where Ai and Bi are suitable constants such that 


Bi = 


Ti'i'iyF) 


Ai. 


Moreover, the eigenvalue A(q associated with the eigenfunction ui is delivered by formula 


\i) = M (1 - yyy/r) + Ti”{y/r 


i{l-l)l{l + N 

-{I- l)y^{N -1 + 2NI + 2l{l - 2)1 + T)i'i{yff) 

+ t{{1 - !)(/ + 2fV - 3) + T)i'{{yy) 

+ {l- l)Ty/Ti'l'{y/r) 


(5.3) 


for any I G N. 

Proof. Solutions of problem (HU in the unit ball are smooth (see e.g., [lU Theorem 2.20]). We 
consider two cases: Au = 0 and Au ^ 0. 

Let u be such that Au = 0. The Laplacian can be written in spherical coordinates as 


A — dyr H- Qt H— 7 A 5 . 

r r^ 

Separating variables so that u = R{r)Y{9) we obtain the equations 


R" 


N-1 , l(l + N-2) 

- R' - ^- ^-R = 0 


(5.4) 


and 

AsY =-l{l + N - 2)Y. (5.5) 

The solutions of equation (15.41) are given by R{r) = ar’" + br^~^~^ if / > 0,iV > 2, and by 
R{r) = a + 61og(r) if Z = 0, A = 2. Since the solutions cannot blow up at r = 0, we must 
impose 6 = 0. The solutions of the second equation are the spherical harmonics of order 1. 
Then u can be written as 

u{r,9) = air'-Yi{9) 

for some I G N. 

Let us consider now the case Au ^ 0. We set v = Au and solve the equation 


Av = TV. 


21 








By writing v = R{r)Y(6) we obtain that R solves the equation 


i?" 




(5.6) 


while Y solves equation ()5.5|) . Equation (jh.BIl is the modified ultraspherical Bessel equation 
that is solved by the modified ultraspherical Bessel functions of first and second kind q(i/rr) 
and ki(y/Tr). Since the solutions cannot blow up at r = 0, we must choose only ii(z) since ki(z) 
has a singularity at z = 0. Then 


v{r,9) = biJi^{^r)Yi^{e) 

for some h € N. Now v = ^ = Au, that is A{v/t — u) = 0. This means that 

u(r,0) = —ii^{y/Tr)Yi^(e) - ciy^Yi^{e) (5.7) 

r 

for some I 2 G N. 

Now we prove that the indexes h and I 2 in (I5.7I) must coincide. This can be shown by 
imposing the boundary condition ^ = 0, which can be written as 

biy[{^)Yye) - ciMh - mM = o. (5.8) 

If the two indexes do not agree, the coefficients oi Yi.^i = 1,2 must vanish since spherical 
harmonics with different indexes are linearly independent on 9f2. Since i" (i/t) > 0, this 
implies = 0 and therefore I 2 = 0 or ^2 = 1- Then we have 

ui{r,9) = (^Ay + Biii{yTr)^Yi{9), (5.9) 

with suitable constants Ai,Bi. In the case ^ 0,1, again from the boundary condition (15.81) we 

have 

lil-l)Ai+Ti''iy^)Bi=0, (5.10) 

then Bi = Ai. Note that the formula holds also in the case / = 0,1 since these indexes 

correspond to i?; = 0. 

Finally, let us consider the boundary condition 



Using in (j5.11|) the representation of ui provided by formula (I5.9|l . we get 


( - A + ;((/- 1)(1 + iV - 2) + t)) Ai + ( - (31(/ + TV - 2) + A)ii(v^) 

- V^{{N -1-2NI- 2{l - 2)1 - T)i'i{V^) + {N- l)^/7i;'(v^) 


+ Ti!!'\yjT 


B, 


Yi{9) = \{Ai + Biii{y^))Yi{9). 


Using equality (15.101) we get that ui given by (15.91) is an eigenfunction of (11.11) on the unit 
ball. Moreover, as a consequence, we also get formula (15.31) for the associated eigenvalue. This 
concludes the proof. □ 


We are ready to state and prove the following theorem concerning the first positive eigen¬ 
value. 

Theorem 5.12. Let LI be the unit ball in centered at the origin. The first positive eigenvalue 
of ([OP is X 2 = A(i) = T. The eorresponding eigenspace is generated by {a:i,a: 2 , ...a:Ar}. 

Proof. By Theorem 15.21 0 = A(o) < r = A(i). We consider formula (15.31) with I = 2. We have 


A(2) = 2(ji'yyf) - 2^2(v^)) dNyVr) - yr{5N -1 + 

+ t{2N -1+ T)i'yy/f) -b Ty/ri'^iy/f) 


(5.13) 
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In order to prove that A( 2 ) > t, we use some well-known recurrence relations between ultras- 
pherical Bessel functions (see [TJ p. 376]), 


i'liVr) = +ii+iiVT), 

(vt) = - -j= -*i(VT)H--Z;+2 (Vt) + *Z-H3(Vt)- 

tsJt t Vt 

Using these relations in (I5.13p . we obtain an equivalent formula for A( 2 ), 

A(2) = 2^5-yTi3(v^) + TH{VT)j (lOA^ -2 + 2t)z2(v^) 

+ {2- fOlV -h (7 -b 107V)^/7 - 2r -b 5rv^)i3(v^) 

-b t( 8 -b 2N -b r)i4(-/r) -b r-\/ri5(-\/r) 


By well-known properties of the functions (see (TJ §9]), it follows that ii > i;+i for all I G N. 
This implies 


(lOiV - 2 -b 2T)i2{VT) + (2 - lOiV -b (7 -b 10fV)VT - 2r -b 5Tv^)i3(v^) 

-b t(8 -b 2iV -b T)i4(i/r) -b r-/ri5(i/r) > (^5r-/ri3(v^) -b T^i4iVT)'j , 


then 

A( 2 ) > 2r > t = A(i). 

Now it remains to prove that A(;) is an increasing function of I for I > 2. We adapt the 
method used in [H Theorem 3]. We claim that for any smooth radial function R{r) the 
Rayleigh quotient 


Q{R{r)Yi{e)) 


\D^iRir)Yi{0))\^ + r|V(R(r)D(g))pdx 
Jg^R{rfYi{eyda 


is an increasing function of I for I > 2. We consider the spherical harmonics to be normal¬ 
ized with respect to the L^{dB) scalar product. In particular, we have that the denominator 
D[R{r)Yi{6)] of Q{R{r)Yi{0)) is R^(I). For the numerator N[R(r)Yi(9)] of the Rayleigh quotient 
we have 


N[R{r)Yi{e)] 



k{k-N-1/2) 


R^ + t- 


kR^ 


„N-1 


dr 



((R"") + ^^{R')^ +T{R')^y^-^dr, 


where k = l{l + N — 2). The above expression is increasing in k for k > N + 1/2 and since k is 
an increasing function of /, we easily get that each term involving I is an increasing function of 
I for I >2. Thus the claim above is proved. 

For each I G N, 


A(;) = inf Q(u) = inf 


-b rjVupda; 

IdB 


(5.14) 


where the infimum is taken among all functions u that are L^(9R)—orthogonal to the first m— 1 
eigenfunctions Ui and m £ N is such that Ap) = A^ is the to— th eigenvalue of problem (11.11) . 
The eigenfunctions ui are of the form ui = Ri{r)Yi{6), and ui realizes the infimum in (|5.14l) . 
Then 


Ap) = Q{Ri{r)Yi{9)) < Q{Ri+4{r)Yi{9)) < Q(R,+i(r)ll+i(0)) = Ap+i), 
where the first inequality follows from the fact that i?/+i(r)Y/(0) is also orthogonal with respect 
to the L^{dB) scalar product to the first to — 1 eigenfunctions Ri{r)Yi{9) for i = 1, ...to — 1, 
and then it is a suitable trial function in (15.141) . The second inequality follows from the fact 
that the quotient Q{R{r)Yi{9)) is an increasing function of I, for I > 2. This concludes the 
proof. □ 
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5.2 The isoperimetric inequality 

In this section we prove the isoperimetric inequality (15.11) . Actually, we prove a stronger result, 
that is a quantitative version of (15.11) . We adapt to our case a result of [8], where the authors 
prove a quantitative version of the Brock-Weinstock inequality for the Steklov Laplacian. We 
also refer to miiia where these kind of questions have been considered for the first time (see 

also mm)- 

Throughout this section 17 is a bounded domain of class . We recall the following lemma 
from [8]. 

Lemma 5.15. Let 17 be an open set with Lipschitz boundary and p > 1. Then 

XoL- (t) )’ 

where 17* is the ball centered at zero with the same measure as 17, 17A17* is the symmetric 
difference of Q and 17*, and cn,p is a constant depending only on N and p given by 


cn,p ■= 


iN+p-l)ip-l) V2-1 


N 


min fP- 


We also recall the following characterization of the inverses of the eigenvalues of (ED from 
pS] (see also [7]). 

Lemma 5.16. Let 17 he a bounded domain of class inM.^. Then the eigenvalues of problem 
EP on ri satisfy, 


fe+iV 

E 


Ai(17) 


= max < 


fc+iV 

E 

Efc + 1 ' 
fc+iV 


an 


vf da }, 


(5.17) 


where the maximum is taken over the families {u/};Xfc-i-i satisfying D^Vi : D^Vj -I- 

rVui • Vvjdx = Sij, and ViUjda = 0 for all i = k + 1, ..., k + N and j = 1, 2,..., k, where 
ui,U 2 , ■■■,Uk are the first k eigenfunctions of problem EP- 

For every open set 17 G with finite measure, we recall the definition of Fraenkel asym¬ 
metry 


A(17) :=inf|E)?—^ ball with |B| = |17|| . 

The quantity A(17) is the distance in the L^(R^) norm of a set 17 from the set of all balls 
of the same measure as 17. This quantity turns out to be a suitable distance between sets 
for the purposes of stability estimates of eigenvalues. Note that A(17) is scaling invariant and 
0 < A(17) < 2. 

We are ready to prove the following 

Theorem 5.18. For every domain 17 in of class the following estimate holds 

A2(17) < A2(17*) (l-<5iv^(17)2) , (5.19) 


where Sn is given by 


Sn ■= 


N + 1 
8N 



and 17* is a ball with the same measure as 17. 


Proof. Let 17 be a bounded domain of class in R'^ with the same measure as the unit 
ball B. We consider in (|5.17l) I = 2,...,N 1 and vi = (t|17|)“^/^ Xi as trial functions. The 

trial functions must have zero integral mean over 917. This can be obtained by a change of 
coordinates x = y — fg^ yda. Moreover, the functions vi satisfy the normalization condition 

of Lemma [5.16l Then vi are suitable trial functions to test in formula ()5.17|) . We get 


Af-l-l y 


> 


-|17| 


l^da. 


'an 
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We use Lemma [5.151 with p = 2. This yields 


. r.+" f 1+c„ r^i 


t\B\ 


I |f^l ) 


+ N J 


AT+l 


= E 


1^2 


MB) 


1 + Cat 9 


f\nmy" 
\ M J ^ 


Suppose now that A 2 (n) > otherwise estimate (15.191) is trivially true, since 0 < Al(n) < 2. 
Since A 2 (n) < A; (11) for all I > 3, the previous inequality and the definition of yl(n) yield 

A2(11) (l + C7V^2Al(n)^) < X2{B). 

This implies (15.191) with 5n = | min{ 1,— !)}■ Note that min{l,^^^^( ^^2 — 1)} = 
^2 — 1). This concludes the proof in the case 12 has the same measure as the unit ball. 
The proof for general finite values of |r2| relies on the well-known scaling properties of the 
eigenvalues. Namely, for all a > 0, if we write an eigenvalue of problem (El as A(r, 12), we 
have 

A(t, 12) = a^A(a“^T, all). 

This is easy to prove by looking at the variational characterization of A(t, 12) and A(a“^r, all) 
and performing a change of variable x x/a in the Rayleigh quotient (13.91) . This last obser¬ 
vation concludes the proof of the theorem. □ 

The isoperimetric inequality 15.11 is an immediate consequence of Theorem 15.181 

Corollary 5.20. Among all bounded domains of class M with fixed measure, the ball maximizes 
the first non-negative eigenvalue of problem EP, that is A 2 ( 11 ) < A 2 (ll*), where A 2 (ll) has been 
defined in iS. 9\) and 12* is a ball with the same measure as 12. 

Remark 5.21. In JE! the authors prove that the quantitative version of the Brock-Weinstock 
inequality that they find is sharp. We think that it would be of interest to consider the problem 
of the sharpness of inequality as well. Unfortunately, the results of JEH do not apply imme¬ 
diately to our case. Also, we do not discuss the sharpness here since we think it is out of the 
purposes of the present paper. Such a discussion will be part of a future work. 


6 Concluding remarks 


Throughout this paper we have only considered problems (II. 1|) and (13.101) with r > 0. If we set 
r = 0, problem (El reads 


( A^u = 0, in 12, 

/ ^ = 0 on (912 

1 -diyon’iDMiy) - ^ = Xu, on 9111 


while problem (13.101) reads 

{ A^u = Xu, in 12, 

0=0, on 911, 

divan -I- = 0, on 912. 


( 6 . 1 ) 


( 6 . 2 ) 


Problems dSH) and (16.21) model free vibrating plates which are not subject to lateral ten¬ 
sion. These problems have a sequence of non-negative eigenvalues of finite multiplicity and the 
corresponding eigenfunctions form a orthonormal basis of iL^(ll). The coordinate functions 
Xi, ...,xn and the constants are eigenfunctions of both problems (16. ip and (16.21) corresponding 
to the eigenvalue A = 0, which has multiplicity N -\-1. Therefore, the first non-zero eigenvalue 
is the {N -\- 2)-th eigenvalue. 
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(6.3) 


As we did in Theorem 13.171 we can define the family of problems 

{ A^it = \pe M, in n, 

0 = 0, on on, 

divan = 0, on 511, 

where is dehned as in (13.1111 . We have the following theorem, whose proof can be easily done 
adapting that of Theorem 13.171 

Theorem 6.4. Let LI be a hounded domain in of class C^. Let p^ be defined as in 
Let \j{pe) be the eigenvalues of problem on LI for all j G N. Let Xj, j € N denote the 

eigenvalues of problem \6.1\) corresponding to the constant surface density -^i-. Then we have 
lime_>o Xj{pe) = Xj for all j G N. 

It is clear that a discussion similar to that of Section |4] can be carried out for problems 
(lOl and (16.21) as well, by means of a change of the projections , tt^ according to the kernel. 
In particular, all the formulas in Section 2] remain true, by setting r = 0. Then we have the 
following 

Theorem 6.5. Let LI be a domain in Let ^ G be sueh that (j>{Ll) is a ball. Let X be an 

eigenvalue of problem i6.1\) (problem 16'. ill] respectively) in 4>(Ll), and let F be the set of j G N 
such that Xj[(j)] = X. Then Ap^g has a critical point at (f on V{V{(j})), for a/Z s = 1,..., \F\. 

Moreover, for problem (16.11) . it is possible to identify the fundamental modes and the fun¬ 
damental tone on the ball. We have the following 

Theorem 6.6. Let LI = B be the unit ball in . The eigenfunctions of problem 16.11] are of 
the form 

ui{r,e) = {Air^ + Bir^+^)Yii9), 
for I G N, where Ai and Bi are suitable constants such that 

Z(Z-l) 

' (Z + 2)(Z + 1)^'- 


The eigenvalues X(y) of problem 16.il] corresponding to the eigenfunctions ui{r,9) are delivered 
by the formula 


\i) - 


HI -1){N + 2Nl +{1- 1)(2 -1 3Z)) 
1-1 2Z 


The first positive eigenvalue is 


Xn+2 — X(2) — 2 ^ 


and the eorresponding eigenfunctions are 

U 2 { r , 9 ) = (6r^ - 1^2(6*)• 


Proof. The proof is similar to that of Theorem 15.21 from which it differs only for the use of 
biharmonic functions on the ball as solutions of the differential equation = 0. For a 
characterization of biharmonic functions on the ball we refer to [21 |3l [38] . □ 


We have an explicit form for the fundamental tone and for the corresponding eigenfunctions 
in the case of the unit ball which suggests how to construct trial functions for the Rayleigh 
quotient of AAf- 1 - 2 . Unfortunately, if we want to use a function of the form R{r)Y2{9) as a 
test function as we did in Theorem 15.181 we must impose that R(r)Y2{9) is othogonal to the 
constants and to the coordinate functions with respect to the Lf{dL>.) scalar product and we 
can no more obtain this just by translating the domain LI. 

We remark that functions of the form R{r)Y2{9) where R(r) = 6r^ — r'^ for r G [0,1] and 
R(r) = 8r — 3 for r > 1 are suitable trial functions for the annuli. Explicit computations show 
that, for example, in dimension 2 or 3 (where the formulas are less involved), the ball is a 
maximizer among radial domains with a fixed measure. 

The results contained in this section suggest that the ball should be a maximizer also for 
problems (lO) and (16.21) . For what concerns problem (16.21) . a characterization of the funda¬ 
mental tone is still unavaiable. A deeper analysis of problems (16.11) and (16.21) will be part of a 
future work. 
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